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ABSTRACT 



CU Virginis is a rapidly rotating Magnetic Chemically Peculiar star with at 
present unique characteristics as radio emitter. The most intriguing one is the 
presence of intense, 100% circularly polarized radiation ascribed to Cyclotron 
Maser. Each time the star rotates, this highly beamed emission points two times 
toward the Earth, like a pulsar. We observed CU Vir in April 2010 with the 
EVLA in two bands centered at 1450 and 1850 MHz. We covered nearly the 
whole rotational period, confirming the presence of the two pulses at a flux den- 
sity up to 20 mJy. Dynamical spectra, obtained with unprecedented spectral 
and temporal sensitivity, allow us to clearly see the different time delays as a 
function of the frequency. We interpret this behaviour as a propagation effect of 
the radiation inside the stellar magnetosphere. The emerging scenario suggests 
interesting similarities with the auroral radio emission from planets, in particular 
with the Auroral Kilometric Radiation (AKR) from Earth, which originates at 
few terrestrial radii above the magnetic poles and was only recently discovered 
to be highly beamed. We conclude that the magnetospheres of CU Vir, Earth 
and other planets, maybe also exoplanets, could have similar geometrical and 
physical characteristics in the regions where the cyclotron maser is generated. In 
addition, the pulses are perfect "markers" of the rotation period. This has given 
us for the first time the possibility to measure with extraordinary accuracy the 
spin down of a star on or near the main sequence. 

Subject headings: stars: chemically peculiar — stars: individual (CU Virginis) — 
stars: magnetic field — polarization — radiation mechanisms: non-thermal — masers 
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Introduction 



Among Magnetic Chemically Peculiar (MCP) stars (see IWolfla (Il983[ ) for an extensive 
description of this class), CU Vir (HD124224) is one of the most intriguing. It is a rapidly 
rotating star, with one of the shortest ro tational period s of this class (Prot ~ 0.52^^). 



According to the Oblique Rotator Model ( iBabcock 



19491 ). it is characterized by a magnetic 



field with a mainly dipolar topology, with a strength of Bp = 3000 ± 200 gauss at the pole, 
whose axis is tilted of an angle /3 = 7 4° ± 3° with respec t to the rotational one, which in 



turn subtends an angle i = 43° ± 7° (ITrigilio et al 



20001 ) to the direction of the Earth. 



The magnetic field is also responsible for the inhomogeneous distribution of the chemical 
elements in the stellar photosphere, giving rise to photometric and spectroscopic variability. 



CUVir shows continuum radio emission (ILeone et al.lll994j . |2004J ). explained as being 



due to gyrosynchrotron process from e 
According the 3D model developed by 



ectrons spira li ng in a corotating magnetosphere. 



Trigilio et al. 



(120041 1 for MCP stars, a radiatively- 



driven wind interacts with the magnetosphere, dividing it into three regions, namely the 
"inner", "middle" and "outer" magnetosphere. In the "inner" equatorial zone, the kinetic 
energy density of the wind particle never exceeds the magnetic one, B'^/Stt > (l/2)pf^, 
generating a dense equatorial belt. In the outer region, the path of the ionized wind from 
the polar cups is first traced by the dipolar field, then, where B'^/Sti < (l/2)pf^, it is mainly 
radial. In the thin middle magnetosphere, particles evaporating from a ring around the 
magnetic poles propagates along the dipolar field lines, breaking them at the Alfven radius 
(-B^/Stt ^ {l/2)pv'^), generating current sheets where particle acceleration occurs. The 
energetic particles spread into the middle magnetosphere, emitting by the gyrosynchrotron 
mechanism. For CUVir, a weak wind with a mass-loss rate of the order of 10^^^ Mq year^^ 
can explain the almost fiat spectra and rotational modulation. 



The distinctive behaviour of CUVir is the presence of pulses of coherent, 100% 
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circularly polarized, highly directive radio emission at 1.4 GHz ( iTrigilio et al.l 120001 ) . Pulses 
are observed two times every stellar rotation, when the magnetic axis of the dipole is 
perpendicular to the line of sight. This pr ocess has been observed for over ten years 



( iTrigilio et al. 



2000 



2008 



Ravi et al. 



2010l ). and it is interpret ed as Electron Cyc l otroii 



2000 



20081 ). 



Maser Emission (ECME) originating above the magnetic pole ( ITrigilio et al.l 
In the framework of the previous model, the relativistic electrons in a converging flux tube 
reflect outwards when they arrive close to the star, some fraction of them falling into the 
stellar photosphere. This causes a loss cone anisotropy which is responsible for the maser 
emission. 

Another important aspect of the two pulses is that, assuming the emission source is 
fixed, they represent an ideal marker of the stellar rotation. An abrupt variat ion of the 



rotational period of CU Vir has been claimed from optical photometric study by 



Pyper et al. 



(119981 ). A c hange of the period from the radio pulses study has been measured with great 



accuracy by 



Trigilio et al 



(|2008[ ). 



In order to study the spectral and temporal behaviour of the pulses, the magnetospheric 
plasma, the spin down of this star, we observed CU Vir with the EVLA with unprecedented 
spectral and temporal sensitivity in two bands around the wavelength of 18 cm. 



2. Observations and Data Reduction 



201l|) in D 



The observations were carried out with the EVLA (iPerley et al. 
configuration, in two days, on April 23rd, from 03:36 to 10:35 UT, and on April 30th 2010, 
form 01:39 to 08:38 UT. The date of the observations have been coordinated in order 
to cover the whole rotational period of CUVir (about 12 hours) with two 7-hours slots, 
with a minimum overlap of rotational phase. We used all the 22 telescopes equipped with 
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Fig. 1. — Dynamical spectra of the circular polarization (Stokes V) flux density of CU Vir; 
the corresponding heliocentric rotational phase is also shown. Left panel) first day: a broad 
peak, labeled "a", centered at 8 UT is visible; Right panel) second day: two peaks, one at 3 
UT, "b", one at 7:30 UT, "c", are visible. Peak "a" and "b" appear at the same rotational 
phase (0.5-0.6); The peak "c" appears at phase 0.95-1. Vertical lines show that the peaks at 
~ 0.55 occur early at low frequency; conversely, the peak at ~ 0.95 occurs early at high 
frequency. 
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the L-band receivers, tuned in two bands centered at 1452 and 1860 MHz respectively, 
each one with a bandwidth of 128 MHz spht into 64 channels, with a frequency resolution 
of 2 MHz, in full polarization mode. The acquisition on CUVir was alternated with 
those on the point-hke source J1354-0206, used as phase cahbrator, with a duty cycle of 
1™30^ — 9™30^ — 1™30^. The amplitude scale was determined by observations of the primary 
calibrator 3C286 (Fi 4 = 14.58, Fi § = 12.9 Jy). Data reduction and editing were performed 
by using the Common Astronomy Software Applications package (casa). Data 
inspection revealed much radio interference in both bands, occurring in certain blocks of 
time and of channels. Bad data have been edited first in the scans of the two cahbrators. 
Bandpass correction has been eventually determined by using the visibilities of 3C286, then 
flux densities of J1354-0206 have been obtained (F1.4 = 774, F1.8 = 854 mJy). Complex 
gains have then been apphed to the data of CUVir which have been edited after cahbration. 
Stokes 1 and V maps of the regions around the phase calibrator and the target have been 
obtained for the two days. 

In order to get dynamical spectra of our target, visibilities of CUVir have been 
imported into the Astronomical Image Processing System (aips), where the task 
DFTPL computes the direct Fourier Transform at a given position in the sky as a function 
of the time. Stokes V flux densities have been computed for all frequencies, with time and 
spectral resolution of 2 MHz and 4min respectively, giving an rms of about 1 mJy, as 
expected. 

No Stokes I or V variability has been found for J1354-0206, neither for the several 
sources close to CUVir. 
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CU Vir - V Polarization - 1450 IVIHz 




Rotational Phase 



Fig. 2. — Frequency average of the polarized emission as a function of the rotational phase 
for the two bands. 

3. The Cyclotron Maser 



Dynamical spectra are shown in 



using the ephemeris by 



Pyper et al. 



Fig, 



H Heliocentric rotational phases, computed by 



(119981 ) ■ indicate that the peak "a" of April 23rd and 
the "b" of April 30th occur at the same phase (0 ~ 0.55), when the star is oriented in space 
in the same way relatively to the Earth. Fig. |2] shows the polarized flux average in the two 
bands as a function of the rotational phase, stressing the coincidence in phase of "a" and 
"b" events and a significant level of variability between the two observations. All the three 
events last approximately 1 hour each; "c" exhibits a single peak of flux levels of about 20 
mJy; "a" and "b" lower and broader emission, with two components in the "b" event. The 
phase difference between "a" (or "b") an d "c" events i s abou t 0.4. This allows to r ecognize 



the same peaks o f emission reported by 



Ravi et al 



Trigilio et al. 



fcoooh 



Trigilio et al. 



(120081 ) and 



(120101 1. The difference with the previous observations is the rotational phase of 
each peak, not the difference of phase between them. The variability of the pulse shape, as 
well as the amplitude, could be due to some instability of the emission region, as proposed 
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by 



Ravi et a. 



flMutel et al 



(20 



0|). Many parameters enter in the growth rate u) of the cyclotron maser 
2007D such as electron density, velocity distribution (/) and its anisotropy 
[df /dvj_), plasma frequency (z/p) and others. Since the growth rate is a non linear function 
of them, and the intensity of an unsaturated maser is / oc exp(a;), even a small variation of 



one of the above parameters can result in a significant variation of the emerging flux. 



3.1. The change of the rotational period 



by 



An abrupt 



Pyper et al. 



period change of 2.18 seconds, occurring probably in 1985, has been claimed 



(119981 ) from the analysis o 



second has been inferred by (iTrigilio et al 



phot ometric data; a new spin down of about 1 



20081 ) from the phase shift of the coherent radio 



emission betwee n 1998 anc 
inclusion of the Ravi et al 



999. Now, with a new set of points after 12 years, and with the 



(120 lOf ) results, the spin down is confirmed and the increment of 
the rotational period AP is determined with great accuracy by fitting the midpoints of the 
two peaks (circles and line in Fig. |3]): AP = 1.12 ± 0.03 s. If we plot the phase of the peaks 
as a function of the epoch (Fig. [3]), we see that they appear with a regularly increasing 
delay, indicating that the rotational period of CUVir has changed. Assuming that the 
change of period occurred close to the observations of 1998, we determine new ephemeris of 
CUVir: 

HJD = 2450966.3601 + 0.52071601^ 

for JD > 2450966, nominally June 1998. It is worthwhile to stress that this method of 
period determination gives a relative error of 7 x 10^^. 
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3.2. The direction of the beams 



After computing the actual rotational phases cj) with the new ephemeris, it is possible 
to get the orientation of the star and its magnetosphere, then the direction of emission with 
respect to the dipole axis. Assuming a n inclination i = 43 ° of the rotational axis and an 



obliquity (3 = 74° of the magnetic axis (ITrigilio et al. 



20001 ). the angle ip between the line of 



sight and the axis of the dipole is given by: 



here 



sin ip = sm/3 sin i cos(0 — 0o) + cos /3 cos 



0.08 is the phas e delay of the magnetic curve with respect to the light curve 



( iBorra fc Landstreet 



19801). The maximum visibility of the North magnetic hemisphere is 



at 



0.58, while a.t (j) = 0.38 and 0.78 the magnetic axis lies in the plane of the sky. 



By looking at the radio emission as a function of the inclination of the magnetic axis, 
we can see how it escapes from the emitting region, i.e. we can study its directivity. The 
direction of the emission beams at 1450 and 1850 MHz relatively to the dipole axis are 
shown in Fig. SI top panel. Beams at higher frequency are emitted at ~ 4°, while at lower 
frequency at ■?/' ~ 0°, as the visibility of the peaks in Fig. [1] indicates. The emitting beams 
subtend an angle Aip < 10° HPBW. 



T he cyc lotron maser hypothesis in CU Vir has already been proposed by iTrigilio et al. 



mod 



20081 ). To be possible, the emission must originate in a region of relatively strong 
magnetic field (z/b >> z/p), and the plasma densi ty in region above the magnetic pole must 



be relatively low. Following the ECME theory (IMelrose fc Dulk 



19821 ). the direction of 



emission forms an angle 6 < 90° (cos 6 ^ v/c, where v is the upward velocity of the emitting 
electrons) to the magnetic field direction with an aperture A^^ ^ 'v/c), forming an hollow 
cone. The frequency of the radiation is u > s ub, where s is the harmonic number. It is 
probably the second harmonic which is the most efficient one, since emission at the first 
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harmonic is subject to gyromagnetic absorption of nearby thermal plasma at s = 2, and for 
a higher harmonic number the maser has less efficient amplification. Assuming s = 2, the 
emission occurs where u = 2 x 2.8B/G MHz. For our two observing frequencies of 1450 
and 1850 MHz, the corresponding magnetic field is respectively 259 and 330 G and, with 
a polar field strength i?poi = 3000 G, the observed emission arises respectively at i? = 2.26 
and 2.09 -R* from the center of the star. 

There are however some problems with the hollow cone beaming model. In fact, if the 
emission arises from a ring around the polar axis and if the emission is in direction 6 < 90°, 
we should observe much broader peaks, since the inclination of the m agnetic field lines wi th 



respect to the line of sight changes in the emitting rings (see Fig. 4 of iTrigilio et al.l ( l2008l )). 
So the open questions are: 

-why do we not observe broad peaks? 

-why does the angle of emission change with the frequency? 



4. "Auroral" emission 



The cyclotron maser of CUVir shares with the radio emission from planets some 
characteristics, suggesting a common scenario. There are five planets with magnetospheres 
in the Solar System (Earth, Jupiter, Saturn, Uranus and Neptune) known to show intense 
radio emission at low frequency (f rom 10 kHz up to few MHz), originating above the high 



magnetic latitude auroral regions (jZarka 



19981 ). This emission is non-thermal, with high 



brightness temperature exceeding lO^^K, 100% circularly polarized, highly beamed. All 
those characteristics indicate that the emission process is cyclotron maser. Except for 
the case of the Jovian decametric emission (DAM) at 10 — 40 MHz, which is due to the 
lo-magnetosphere interaction, the auroral emissions are controlled by the solar wind fiowing 
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into the open field lines above the polar caps. The analogy with CUVir is unexpected, 
including the "pulsar like" behaviour due to the rotation of an oblique dipole axis. 

Recent observations of the terrestrial Auroral Kilometric Radiation (AKR), performed 
with the four-spacecraft Clust er array, have allo wed the location of the sources from 



which the emission origi nates (iMute 



beaming characteristics ( iMutel et al 



et al. 



20031 ) and the determination of its the angular 



20081 ). The emission seems to be the sum of several 
elementary bursts of radiation generated at a height where the gyrofrequency of the 
electrons corresponds to the magnetic field (at 0.5 — 2R@), in a ring around the magnetic 
axis where the density is very low (the so-called "auroral cavity"). Unlike the Jovian 
DAM, wher e the emission is beamed in a thin (A^^ ^ 1°) hollow cone of half-angle 9 < 90° 



(iDulk 



19701 ). the AKR is radiated within ±15° from the "plane tangent to the source's 



magnetic latitude circle an d containing the local magnetic field vector" (IMutel et al. 



2008 



Louarn &: Le Queau 



19961 ). The important point we want to stress is that the emission is 



emitted tangentially to the auroral ring. 

During the propagation in the close ambient medium, the AKR is deviated upwards 
because of the refractive index of the dense magnetospheric plasma. 



4.1. A new model for CU Vir 



The results of AKR suggest that the ECME of CUVir could come from the same 
geometry. We propose a new scenario, where the masing emission does not follow the 
"hollow cone model", instead the "tangent plane beaming model". Fig. H] shows a sketch of 
the Cyclotron Maser emission from CUVir. Electrons travelin g back from the acceleration 



region located at Alfven radius -Raw = 15 (ILeto et al. 



20061 ) emit at locations where the 
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magnetic field strength corresponds to the first or second gyromagnetic harmonic [s = 1, 2). 
These define annular rings above the magnetic pole. The amplification occurs tangentially, 
so that it is seen from Earth to came from the two extreme points of the rings. They would 
be seen simultaneously at any frequency as the star rotates, in the absence of subsequent 
refraction. 



Our 3-D model (ITrigilio et al.l 120041 : iLeto et al.ll2006l ) foresees the existence of a cold 



torus around the magnetic equator. The absorption of this optically thick material can 
account for the rotational modulation of the continuous gyrosynchrotron emission, in 
particularly at high frequency radiation {u > 15 GHz), emitted exactly in the same region 
where the cyclotron maser originates. From the results of the simulation, the number 
density of the cold torus is iVtorus ~ lO^cm"^, with a temperature Ttorus ~ lO'^K. 

We now consider the case of s = 2, and will give some consideration for the case s = 1. 
During the propagation, the maser radiation crosses the torus. The plasma of the torus, 
permeated by a magnetic field B ^ 260 and 330 G (the levels of the propagation of the two 
frequency bands) has a refractive index in the extraordinary mode given by: 



4 



In the following, we will give an order of magnitude estimate of the deviation of the emission 
(angle ip) with respect to the original, horizontal, ray path. The refractive index inside the 
torus, assuming an average A^toms as above, but with higher density in the inner part, and 
the local field strength, can easily reach values of nrefr ~ 0.98 and 0.95 for u = 1450 and 
1850 MHz respectively. This means that the maximum deviation toward high values of the 
angle ip occurs at 1850 MHz, in agreement with our observations. The actual angle ip can 
be obtained in a simply way from the Snell law if the angle of incidence with the interface 
between the two media at different refractive index, i.e. the surface defined by the field lines 
confining the cold torus, is known. In the case that this angle is 60°, ip ~ 10° at 1850 MHz 



- 13 - 



and ip ^ 5° ai 1450 MHz, accounting for the observed difference (see Fig. H]). 

In the case s = 1, this simple computation foresees bigger deviations upwards, not 
in agreement with the observations. However, a more detailed numerical computation is 
needed in order to better constrain the actual harmonic number s and get more precise 
indications on the physical conditions of the emitting plasma and of the surrounding region. 

5. Conclusion and future perspective 

The observations reported here of CUVir confirm the presence of a steady cyclotron 
maser operating over more than 10 years in the magnetosphere of this star. Our model 
proposes that the process is triggered by the continuous injection of high energy electrons 
flowing from the Alfven region, about 15 stellar radii from the star, to the inner parts of the 
magnetosphere. Those electrons are also responsible for the gyrosynchrotron emission. The 
overall mechanism is powered by the radiation driven wind interacting with the magnetic 
field. 

The high temporal and spectral resolution of the current EVLA allowed us to get 
dynamical spectra over two broad spectral regions. The observations have been carried 
in two days, covering almost the whole rotational period of CUVir. This allowed us to 
measure with great accuracy the delay of the peaks of emission, and so to determine with 
unprecedented precision 10^^) the rotational period of the star, confirming the spin 
down of this young star. Further observations are needed in order to monitor the rotation of 
CUVir, since the emission peaks are believed to be "perfect markers" of this stellar clock. 

There is an impressive similarity between the cyclotron maser of CUVir and the 
coherent radio emission from planets. The recent observations of the Auroral Kilometric 
Radiation from the terrestrial magnetosphere and the new interpretation on the basis of the 



- 14 - 



"tangent plane beaming model" gave the opportunity to formulate a new model for CU Vir, 
where the emission is " auroral- type" , emitted tangentially to the auroral circles. 

The analysis of the dynamical spectra permitted us to clearly see that the pulses 
occurred at slightly different times inside the frequency band. In the framework of the new 
model, this could be due to refractive effects due to the propagation of the radiation inside 
a denser magnetized plasma surrounding the star in the magnetic equatorial belt. This is a 
further indication of the existence of a "cold torus", which also acts as an absorber of the 
gyrosynchrotron continuous emission at centimeter wavelengths. 

In the near future the EVLA will increase its maximum bandwidth up to 2 GHz. 
This will allow us to extend the dynamical spectra in a broader spectral region and to 
better define the low and high limits of the cyclotron maser. The time visibility of the 
peaks as a function of the frequency will allow us to study and model with great precision 
the propagation of the radiation in the magnetospheric plasma surrounding the star. 
CUVir is an excellent astrophysical laboratory for plasma physics in stellar and planetary 
environment. 

We thank the NRAO scheduler who coordinated the times of the two observations to 
cover almost the whole rotational period of CU Vir. The EVLA is a facility of the National 
Radio Astronomy Observatory which is operated by Associated Universities, Inc. under 
cooperative agreement with the National Science Foundation. 

Facilities: EVLA. 
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Fig. 3. — Phases of the peaks at the different epochs of observations, as a function of the 
number of rotations of the star from the first observation. Years are also shown at the top. 
The blue diamonds refer to the main emission beams, red dots to the midpoints. From the 
slope of the linear fit the increment of the rotational period can be inferred. 
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Fig. 4. — Top panel: polar diagrams of the pulses, showing their directivity. At 1450 MHz the 
emission is almost perpendicularly to the magnetic axis B; at 1850 MHz there is a deviation 
of about 4°. Bottom panel: picture of the proposed model: maser amplification occurs in 
annular rings around the pole, at height determined by the gyrofrequency of the emitting 
electrons; the emission occurs tangentially, then refraction in the higher density region (the 
"cold torus" in the inner magnetosphere) deviates the radiation upward; this effect is bigger 
at higher frequency. During the rotation, the pulses are observed from the Earth (supposed 
to be toward the right) two times, at slightly different moments at different frequencies. 
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